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SUMMARY 

Turbine blades for reusable space propulsion systems are subject to severe 
thermomechanical loading cycles that result In large Inelastic strains and 
very short lives. These components require the use of anisotropic high- 
temperature alloys to meet the safety and durability requirements of such 
systems. To assess the effects on blade life of material anlsotrophy, cyclic 
structural analyses are being performed for the first stage high-pressure fuel 
turbopump blade (HPFTB) of the space shuttle main engine (SSME). The blade 
alloy Is directionally solidified MAR-M 246 alloy. The analyses are based on a 
typical test stand engine cycle. Stress-strain histories at the airfoil 
critical location are computed using the MARC nonlinear finite-element 
computer code. The MARC solutions are compared to cyclic response predictions 
from a simplified structural analysis procedure developed at the NASA Lewis 
Research Center. 


INTRODUCTION 

Fuel turbopumps for reusable space propulsion systems operate under 
extreme gas pressure and temperature environments. These operating conditions 
subject the high-pressure stage turbine nozzles and blades to severe thermal 
transients that result In large Inelastic strains and rapid crack Initiation. 
To attain even the short lives required for these systems necessitates the use 
of anisotropic turbine blading alloys. Assessing or Improving the durability 
of these components Is contingent on accurate knowledge of the stress-strain 
history at the critical location for crack Initiation. 

Nonlinear finite-element analysis techniques have become available In 
recent years for calculating Inelastic structural response under cyclic 
loading. These methods are based on classical Incremental plasticity theory 
with uncoupled creep constitutive models. Many of the nonlinear 
finite-element computer codes such as MARC (ref. 1) have the capability of 
handling materials with anisotropic properties. However, these codes are 
usually too costly and time consuming to use In the early design stages for 
aerospace applications. Costs are further Increased by the geometrical 
complexity of high-pressure turbine blades which require three-dimensional 
analyses and sometimes substructuring to obtain accurate solutions. To 
Improve the design of engine hot path components such as turbine blades, 
simplified procedures for more economically representing structural response 
unde;- cyclic loading have been under development (refs. 2 to 4). 

This paper addresses the problems of calculating the structural response 
of high-temperature space propulslor. components such as fuel turbopump turbine 
blades. The first high-pressure stage fuel turbine blade (HPFTB) In the 


liquid hydrogen turbopump of the space shuttle main engine (SSME) was selected 
for this study. In the past these blades have undergone cracking In the blade 
shank region and at the airfoil leading edge adjacent to the platform. To 
achieve the necessary durability, these blades are currently being cast using 
directional solidification. Single crystal alloys are also under Investigation 
for future SSME applications. The objective of the study was to evaluate the 
utility of advanced structural analysis methods In assessing the low-cycle 
fatigue lives of these anisotropic components. MARC elastic-plastic 
finite-element analyses are being performed for the blade airfoil. Because of 
the extensive computation time required for the nonlinear finite-element 
analyses, neither the blade platform nor shank regions were modeled at this 
time. Analyses of the shank region are also hindered by lack of knowledge as 
to Its thermal environment. A simplified Inelastic structural analysis 
procedure will also be exercised on this problem to assess Its applicability 
to anisotropic turbine blades for space propulsion systems. 


PROBLEM 0ESCRIP1 ION 

The turbine blade airfoil of the high pressure stage of the SSME fuel 
turbopump was analyzed because of Its history of rapid crack Initiation. 

These uncooled airfoils have a span length of 2.2 cm and a span-to-chord width 
aspect ratio of approximately unity. The blades are cast from directionally 
solidified MAR-M 246 alloy. Temperature-dependent properties for this alloy 
were mainly provided by the Rocketdyne Division of Rockwell International 
Corporation. Material elastic properties are summarized In table I. Mean 
thermal coefficient of expansion data were converted to Instantaneous values 
for MARC Input. Poisson's ratio values of 0.143 and 0.391 were used for the 
longitudinal and transverse directions, respectively, and were assumed to be 
constant with temperature. Longitudinal stress-strain properties, summarized 
In table II, were used for the elastic-plastic region. A single crystal MAR-M 
246 alloy Is also being considered for turbine blades In future SSME 
applications. 

Cracking has occurred during service at the airfoil base near the leading 
edge and In the blade root shank area. These cracks apparently Initiate 
during the first few mission cycles due to the severe thermal transients and 
are propagated by vibratory excitation. Since the primary purpose of this 
study was to compare nonlinear finite-element and simplified analytical 
methods, the blade root and platform were excluded rrom the analysis to reduce 
the size of the problem and, therefore, the computing time. 

The mission used for the analysis Is shown In figure 1 In terms of turbine 
Inlet temperature, gas pressure and rpm. This cycle Is applicable to a 
factory test of the engine; It Is also reasonably representative of a flight 
mission except for the foreshortened steady-slate operating time. The major 
factor Inducing rapid cracking Is the transient thermal stresses and the 
Inelastic strains caused by the sharp acceleration and deceleration transients. 

Rudimentary transient and steady-state three-dimensional heat transfer 
analyses have been conducted using the MARC code. Film coefficients were 
obtained from preliminary Information supplied by Rocketdyne. The gas sink 
temperature was assumed constant around the airfoil surface for each time 
step. Colder boundary conditions were assumed at the airfoil base to simulate 
the effects of the cooling of the blade-to-dlsk attachment region by the 
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liquid hydrogen fuel. Improved transient heat transfer analyses using more 
exact boundary conditions will be performed In the near future to obtain more 
accurate blade temperature Input for the nonlinear structural analysis. 


ANALYTICAL PROCEDURE 

Elastic-plastic analyses are being conducted for the HPFTB airfoil with 
both the MARC code and a simplified analytical procedure developed at the NASA 
Lewis Research Center. The MARC analysis used the anisotropic material 
properties while the simplified analysis used only the Young's modulus and 
Poisson's ratio for the longitudinal (spanwlse) direction. 


MARC Finite-Element Analysis 

A three-dimensional finite-element model of the airfoil (fig. 2) was 
constructed of eight-node Isoparametric elements. The model consisted of 360 
elements with 576 nodes and 1661 unsuppressed degrees of freedom. This model 
was a shortened version of a finite-element model created by Lockheed for a 
NASTRAN elastic analysis of the HPFTB blade at steady state condition; this 
work was done under contract to the NASA Marshall Space Flight Center. The 
main difference between the NASA Lewis and Lockheed models was that the blade 
base and most of the platform were omitted for the MARC nonlinear analysis to 
reduce the computing time and to run the problem In-core on the CRAY computer 
system at Lewis. Boundary conditions were applied to constrain all nodes at 
the base of u.° i.iodel to lie on a plane. Additional boundary conditions were 
Imposed to prevent rigid body motion In this plane. 

The MARC code has been used extensively at Lewis for Inelastic analyses of 
aircraft turbine blades and combustor liners and of space power components. 

In conducting a cyclic analysis, the loedlng history Is divided Into a series 
of Incremental load steps which are sequentially analyzed. The plasticity 
algorithm Is based on a tangent stiffness approach In which the stiffness 
matrix Is reformulated and reassembled for every plastic load Increment. The 
Incremental loads are modified by residual load correction vectors to Insure 
that the solution does not drift from a state of equilibrium. Convergence for 
the Iterative plasticity analysis Is Indicated when the strain energy used In 
assembling the stiffness matrix approximately equals the energy change 
resulting from the Incremental solution. 

The mission cycle was subdivided Into 129 Increments. Incremental loading 
Included centrifugal and gas pressure loads and metal temperature 
distributions as calculated from the heat transfer analysis. Approximately 
1000 words of core storage on the CRAY system were needed to run the problem. 
Each cycle of analysis required about 15 hr of computing time on the Cray 
system. 

The directionality of the elastic material properties causes anisotropic 
constraints. Lekhnltskll (ref. 6) has derived the generalized elastic strain 
equations for an anisotropic body with a transverse plane of Isotropy. Matrix 
Inversion of these equations to solve for the stresses results In the 
relationship 
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where n = E'/E and a = nE/((l ♦ w)(l - v - 2n«' 2 )). Here E ' , G' t and 
denote the Young's moduVi' shear modulus and Poisson's ratio, respectively, 
for the longitudinal or span direction while E, G, and v denote these 
constants with respect to any direction In the transverse plane of Isotropy. 
This anisotropic stress-strain law was Incorporated In the MARC user 
subroutine, HOOKLW. 

Plastic strain calculations were based on Incremental plasticity theory 
using the von Mlses yield criterion, the normality flow rule and a kinematic 
hardening model. The material elastic-plastic behavior was specified by the 
yield strengths and work hardening properties In the longitudinal direction; 
transverse properties were not available. Creep analyses are not being 
performed at the present time because of Inadequate knowledge of the creep 
characteristics for the anisotropic blade material. 


Simplified Analysis 

The simplified analytical procedure was developed to economically 
calculate the stress-strain history at the critical fatigue location of a 
structure subjected to cyclic thermomechanical loading. This procedure has 
been exercised on a wide variety of problems Including multlaxlal loading, 
nonlsothermal conditions, different materials and constitutive models, and 
dwell times at various points In the cycles. Comparisons of the results of 
the simplified analyses with MARC Inelastic solutions for thesp problems have 
shown reasonably good agreement. 

The basic assumption Is that the total strain ranges calculated from 
linear elastic and nonlinear analyses are approximately equal and, therefore, 
the material cyclic response can be calculated using as Input the total strain 
history obtained from an elastic analysis. This assumption Is essentially 
true for thermally dominated loading. A corollary of this assumption Is that 
the stress-strain hysteresis loop calculated from an elastic- plastic analysis 
will be of similar size and shape to that calculated from an elastic analysis, 
but translated In the stress dlrectlor. There Is a version (ref. 4) of the 
procedure that uses Neuber-type corrections (ref. 7) to account for strain 
redistribution due to mechanical loading; however, this version was not 
utilized for this study because of the dominance of the thermal loading. 
Classical Incremental plasticity methods are used to characterize the yield 
surface by a yield condition to describe yielding under multlaxlal stress 
states and by a hardening model to establish the location of the yield surface 
during cycling. This procedure can accommodate Itself to any yield criterion 
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or hardening model. The only requirements are that the elastic Input data, 
whether calculated or measured, be In a form consistent with the yield 
criterion and that the appropriate material properties be used In conjunction 
with the hardening model. 

Ab in most nonlinear computer codes, the von Mlses yield criterion has 
been used In applying stress-strain results from elastic finite- element 
analyses of muitlaxlal problems as Input for the simplified procedure. To 
compute cyclic hysteresis loops for life prediction purposes, the Input von 
Mlses stresses and strains have to be assigned signs, usually on the basis of 
the signs of the dominant principal stresses and strains. 

The elastic Input data are subdivided Into a sufficient number of 
Increments to define the stress-strain cycle. These Increments are analyzed 
sequentially to obtain the cumulative plastic and creep strains and to track 
the yield surface. An Iterative procedure Is used to calculate the yield 
stresses for Increments undergoing plastic straining. First, an estimated 
plastic strain Is assumed for calcclailng an Initial yield stress from the 
stress-strain properties and the simulated hardening model. Then a new 
plastic strain Is calculated as the difference between the total anu elastic 
and creep strain components. The yield stress Is then recalculated using the 
new plastic strain. This Iterative procedure Is repeated until the new and 
previous plastic strains agree within a tolerance of 1 percent. Creep 
computations are performed for Increments Involving dwell times using the 
creep characteristics Incorporated In the code. Depending on the nature of 
the problem, the creep effects are determined on the basis of one of three 
options to be selected; (1) stress relaxation at constant strain, (2) 
cumulative creep at constant stress, or (3) a combination of (1) and (2). 

A FORTRAN IV computer program ( ANSYMP) was created to automatically 
Implement the simplified analytical procedure. A detailed description of the 
calculatlonal scheme Is presented In previous papers on the development of 
this procedure. 


DISCUSSION OF RESUL1S 

Calculated metal temperatures at the leading edge at midspan and at the 
crack Initiation site at the base of the airfoil (critical location) are 
presented In figure 3 as a function of elapsed time during the cycle. The 
assumed gas sink temperature around r.he airfoil Is also plotted. Of 
particular note Is that the leading edge temperature at the airfoil base Is 
cooler than at midspan throughout the cycle. This seems reasonable because of 
the cooling of the blade-to-dlsk attachment region by the liquid hydrogen 
fuel. The colder airfoil base temperatures Induce tensile thermal stresses at 
the critical leading edge location that are additive to the centrifugal 
stresses . 

At the beginning of this study, a steady-state MARC elastic-plastic 
analysis of the airfoil was conducted for rated power level (RPL). The 
temperature Input was obtained from a NASTRAN steady-state heat transfer 
analysis performed by Lockheed. In contrast to the large temperature 
gradients calculated from the NASA heat transfer analysis, Lockheed had a 
maximum temperature variation In the airfoil of only 35 °C. Effective stress 
distributions as calculated at the Gaussian Integration points closest to the 
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suction and pressure surfaces are presented In figure 4. Since the airfoil 
was essentially at an Isothermal state, the stresses primarily reflect the 
centrifugal and gas pressure loading and were entirely In the elastic range. 
The maximum effective stress was only 40 percent of the yield strength. It Is 
apparent that the analysis would not show a low-cycle fatigue problem If 
transient thermal effects were not considered. A MARC nonlinear finite- 
element analysis for the complete mission cycle has been Initiated but, as of 
this writing, has not been completed. 


CONCLUDING REMARKS 

This work Is still In a preliminary stage. A tran: ent heat transfer 
analysis of the HPFTB airfoil has been completed. This analysis Is 
rudimentary because of the Inadequate definition of the heat transfer 
coefficients and omission of the cooled shank region In the finite-element 
model. However, It Is sufficient for the purpose of exercising and evaluating 
the structural analysis methods. A MARC cyclic stress-strain analysis has 
been Initiated. 

The simplified procedure will probably require further development for 
application to SSME turbopump blade problems. This Is not only because of the 
use of anisotropic materials, but also because of the nonproportlonallty of 
the thermal and mechanical loading during the cycle. Experience has Indicted 
that effective stresses and strains are not suitable criteria for relating 
multlaxlal nonproportional stress-strain states to uniaxial material 
properties. 
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TABLE I. - OS MAR-M 246 PHYSICAL PROPERTIES 


Temperature, 

°C 

Modulus of elasticity 
GPa 

Means coefficient of 
thermal expansion, 
percent/°C 

Longitudinal 

Transverse 

21 

131 

183 


93 

120 

179 

0.00113 

204 

125 

175 

.00130 

316 

124 

173 

.00133 

427 

119 

166 

.00141 

538 

114 

162 

.00148 

649 

109 

156 

.00149 

760 

103 

149 

.00156 

671 

97 

142 

.00160 



TABLE II. - DS MAR M 246 STRESS- 
STRAIN PROPERTIES 
(LONGITUDINAL) 







(a) Turbine inlet temperature. 



(b) Turbine inlet pressure. 
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(c) Blade rotational speed. 

Figure 1. - Mission cycle used for analysis. 
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Figure 2. - Airfoil finite element 
model. 
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Figure 3. - Airfoil temperature cycle. 




